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Quantitative study of the spin Hall magnetoresistance in ferromagnetic 
insulator/normal metal hybrids 

Matthias AlthammerJ'^ElSibylle Meyer,'' Hiroyasu Nakayama,'^'"* Michael Schreier,'' Stephan 
Altmannshofer,^ Mathias Weiler,^ Hans Huebl,^ Stephan Geprags,'' Matthias Opel, ^ Rudolf Gross, ^'^ 
Daniel Meier, ^ Christoph Klewe,^ Timo Kuschel,^ Jan-Michael Schmalhorst,^ Gunter Reiss,^ Liming 
Shen,^ Arunava Gupta, ^ Yan-Ting Ghen,'^ Gerrit E. W. Bauer, ^■^''^ Eiji Saitoh, ^'^■■'°'^ and Sebastian T. B. 
Goennenwein'''0 

'Walther-Meissner-Institut, Bayerische Akademie der Wissenschaften, Walther-Meissner-Strasse 8, 85748 Garching, 
Germany 

'University of Alabama, Center for Materials for Information Technology MINT and Dept. Clieni, Tuscaloosa, 
AL 35487 USA 
~ ' Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 

' Laboratory for Nanoelectronics and Spintronics, Research Institute of Electrical Communication, Tohoku University, 
Sendai 980-8577, Japan 

' Physik-Department, Technische Universitdt Miinchen, 85748 Garching, Germany 
^^ Fakultdt fiir Physik, Universitdt Bielefeld, 33615 Bielefeld, Germany 

''^Kavli Institute of NanoScience, Delft University of Technology, 2628 C J Delft, The Netherlands 
^WPI Advanced Institute for Materials Research, Tohoku University, Sendai 980-8577, 
Japan 

'CREST Japan Science and Technology Agency, Tokyo 102-0076, Japan 
'The Advanced Science Research Center, Japan Atomic Energy Agency, Tokai 319-1195, 
Japan 

(Dated: 24 April 2013) 

We experimentally investigate and quantitatively analyze the spin Hall magnetoresistance effect in ferromagnetic insu- 
lator/platinum and ferromagnetic insulator/nonferromagnetic metal/platinum hybrid structures. For the ferromagnetic 
insulator we use either yttrium iron garnet, nickel ferrite or magnetite and for the nonferromagnet copper or gold. 
The spin Hall magnetoresistance effect is theoretically ascribed to the combined action of spin Hall and inverse spin 
Hall effect in the platinum metal top layer It therefore should characteristically depend upon the orientation of the 
magnetization in the adjacent ferromagnet, and prevail even if an additional, nonferromagnetic metal layer is inserted 
between Pt and the ferromagnet. Our experimental data corroborate these theoretical conjectures. Using the spin Hall 
magnetoresistance theory to analyze our data, we extract the spin Hall angle and the spin diffusion length in platinum. 
For a spin mixing conductance of 4 x lO'^ il^'m^^ we obtain a spin Hall angle of 0.11 ± 0.08 and a spin diffusion 
length of ( 1 .5 ± 0.5) nm for Pt in our thin film samples. 



PACS numbers: 72.25.Mk,72.25.Ba,75.47.-m 



I. INTRODUCTION 

Pure spin currents - which transport only (spin) angu- 
lar momentum, but no electrical charge - represent a new 
paradigm for spin transport and spin electronics. In the 
last few years, two spin current generation methods have 
evolved: the spin Seebeck effect (SSE)i"— and the spin 
pumping effecti^"— . Both methods involve ferromagnet 
(FM)/nonferromagnet (NM) hybrid structures to generate and 
to detect the spin currents. While originally the focus was on 
electrically conducting FM/NM heterostructures, ferromag- 
netic insulators (FMI) are increasingly exploited in FMI/NM 
hybrids] ^"i^^ FMIs are electrically insulating materials which 
exhibit long-range magnetic order, such that magnetic exci- 
tations (spin currents) can propagate in FMIs, while charge 
currents cannot. In this sense, FMIs allow to cleanly sepa- 



rate spin current from charge current effects. One of the pro- 
totype examples for a FMI compound is yttrium iron garnet 
(Y3Fe50i2, YIG),i2i2128 

The interplay between spin and charge transport in FM/NM 
devices gives rise to interesting physical phenomena. A 
prominent example is the spin Hall magnetoresistance (SMR) 
discovered recently in FMI/NM hybrids.^"ii The SMR is re- 
lated to the absorption/reflection of a spin current density Js 
flowing along the direction normal to the FMI/NM interface. 
The spin current is generated by a charge current density Jq in 
the NM layer via the spin Hall effect (SHE)i^2i^ 
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The total Js is here the direction of the spin current, while 
its spin polarization s is oriented perpendicular to Js and Jq. 
0!sH = Osh/c is the spin Hall angle defined by the ratio of 
the spin Hall conductivity cJsh and the electric conductivity 
C7 1 34] and e is the positive elementary charge. In FMI/NM 
hybrids, the amount of spin current absorption/reflection by 
the FMI at the FMI/NM interface depends on the orientation 



of the FMI magnetization M with respect to the polarization 
s of the spin current. The amount of spin current reflected 
at the interface in turn induces a charge current via the in- 
verse spin Hall effect (ISHE) in the NM layer The total M- 
orientation dependent spin current flow across the FMI/NM 
interface represents a dissipation channel for charge transport, 
and thus affects the resistance of the NM. The resulting spin 
Hall magnetoresistance of the NM is clearly discernible from 
a conventional anisotropic magnetoresistance-^ (AMR) type 
of effect, since the SMR depends on the angle between the 
magnetization M and the spin polarization s _L Jq, while the 
AMR depends on the angle between M and the charge current 
direction Jq. 

Here, we give a detailed, quantitative analysis of the SMR 
in YIG/Pt, YIG/Au/Pt, and YIG/Cu/Pt samples. We fur- 
thermore report the observation of SMR in nickel ferrite 
(NiFe204)/Pt, and magnetite (Fe304)/Pt samples. The latter 
two materials are also important for SSE measurements since 
recently the SSE was observed in semiconducting NiFe20/f^ 
and weakly conducting Fe30it^. 

The paper is organized as follows: First, we explain the 
SMR effect in a simple phenomenological picture, followed 
by a quantitative model. Next, we present structural and mag- 
netic data of the YIG/Pt hybrids, in connection with exten- 
sive magnetoresistance measurement data as a function of the 
magnetization orientation m. We use our theoretical model 
to quantitatively analyze the experimental data and thus ob- 
tain estimates for the unknown parameters, i.e., the spin Hall 
angle agn and spin diffusion length Apt in Pt. Moreover, we 
verify that the SMR does not originate from a static equilib- 
rium magnetic proximity effect in the non-magnetic Pt layer, 
since the SMR persists even when a second, diamagnetic NM 
is inserted between YIG and Pt. Finally, we show that the 
SMR effect is also observed in FMI/Pt hybrid structures using 
Fe304 or NiFe204 as the FMI. 



II. THEORY OF THE SPIN HALL MAGNETORESISTANCE 

The SMR effect stems from a combination of the spin Hall 
effect (SHE), which converts a charge current Jq into a spin 
current Js, and the inverse spin Hall (ISHE) effect, where a 
spin current Js induces a charge current Jq .^^'^^ 

We start with a phenomenological explanation of the SMR 
effect by discussing the influence of two different boundary 
conditions on the steady state in case of both the ordinary Hall 
effect (OHE) and the SHEi^ We first consider the OHE in the 
single band model with a longitudinal charge current density 
Jq in —X direction as illustrated in Fig.dta). Due to an external 
magnetic field B applied along z, the charge carriers experi- 
ence a Lorentz force in the transverse (y) direction. In typical 
Hall effect measurements, the potential drop (voltage) along y 
is recorded at zero current flow along y. This corresponds to 
the boundary condition that the charge carriers can not escape 
and therefore accumulate on the transverse sides (open cir- 
cuit condition for charge transport). The resulting charge ac- 
cumulation leads to a compensating electrical field Enaii || y, 
which can be detected as a voltage drop over the transverse 




FIG. 1 . (Color online) (a) Illustration of the ordinary Hall effect in 
the single band model. The applied external magnetic field B || z 
deflects the longitudinal charge current Jq to one side of the sample 
due to the Lorentz force, (b) An open circuit boundary condition for 
charge transport on the transverse sides leads to a charge accumu- 
lation on the transverse sides of the sample, generating a transverse 
electric field Ey. It compensates the Lorentz force and leads to a 
steady state in which the longitudinal resistance does not depend on 
B. (c) Short-circuiting the transverse sides results in a transverse 
current flow Jq.y. In steady state, the longitudinal resistance then 
depends on the applied external magnetic field, (d) Due to the spin 
Hall effect a charge current Jq induces a spin current Js || — y with 
a spin polarization s || z (spin up and spin down charge carriers are 
deflected in opposite directions), (e) For open circuit conditions for 
spin transport on the transverse sides (no transverse flow of spin cur- 
rent) a gradient in the spin dependent electrochemical potential fl^ 
(spin accumulation) is generated, such that the spin diffusion current 
compensates the spin Hall currents. In this steady state the longitudi- 
nal charge current resistance does not depend on the spin Hall effect, 
(f) A transverse spin current short circuit (a transverse spin potential 
short), in contrast, results in an additional spin current in the trans- 
verse direction. This effectively leads to a change in the longitudinal 
electrical resistance, due to the inverse spin Hall effect. The change 
in longitudinal electrical resistance between (e) and (f) is the SMR 
effect. 



sides. In the single band model, one finds that the longitudi- 
nal resistance is independent of B (cf. Fig. [Ub))!^ In con- 
trast, one may assume the boundary condition that Enaii = 
(placing an electrical short between the transverse contacts 
as illustrated in Fig. [TJc)), such that a transverse current Jq.y 
arises, depending on the applied magnetic field magnitude B. 
This effectively reduces Jq 1 1 — x. The external current source 
has to compensate for this effect by increasing the effective 



longitudinal voltage Viong applied to the sample, such that 
the charge current /q and thus Jq || — x stays constant. For 
this boundary condition, the observed longitudinal resistance 
^long = Viong/Zq thus does depend on the external magnetic 
field.^ 

We now turn to the SHE illustrated in Fig. [TJd), by which 
the longitudinal charge current density Jq induces a trans- 
verse spin current Js. Since we are here concerned with the 
spin current flow across the FMI/NM interface, we can re- 
strict the discussion to — Js along the interface normal y. Ac- 
cording to Eq.dTJ, this implies that the spin orientation s||z. 
The vector Js thus represents a flow of z-polarized spins along 
— y. For an open (spin current) circuit boundary condition 
as illustrated in Fig. [TJe), we obtain a spin accumulation on 
the transverse sides of the samplers This spin accumula- 
tion leads to a gradient in the spin dependent electrochemical 
potential /is, which compensates the spin current generated 
by the SHE. In the steady state, the longitudinal resistance 
is thus independent of the SHE, in analogy to the OHE in 
Fig. [lib). In contrast, if we consider short-circuited boundary 
conditions for the spin channel (shorting Vjis) a spin current 
flows through the lateral faces. The current source driving 
the charge current again must compensate for this transverse 
current flow, which results in an increase in the longitudinal 
resistance. This change in longitudinal resistance is the SMR 
effect. Dyakonov predicted a magnetoresistance, stemming 
from the magnetic field-dependent dephasing of a spin accu- 
mulation generated by the spin Hall effect at the edges of sin- 
gle NM films4i The spin accumulation at the edges and the 
corresponding magnetoresistance is negligibly small for ex- 
tended thin films. In contrast, we consider here the spin accu- 
mulations on the surface of the conducting film and modulate 
them by a FMI rather than an appHed magnetic field. 

Interestingly, FMI/NM heterostructures allow to gradually 
switch between open and closed spin current circuit boundary 
conditions by changing the magnetization orientation in the 
FMI with respect to s. To realize an ideal spin current short 
circuit condition as depicted in Fig. ^f), FMI/NM/FMI mul- 
tilayer structures are necessary. However, the SMR effect also 
occurs in "simple" FMI/NM hybrid structures, if the spin flip 
length in the NM is finite. We consider first a charge current 
density Jq flowing along j (Fig. |2c)) through the NM layer 
as indicated by the blue arrow in Fig. EJa). This charge cur- 
rent induces a spin current Js || n due to the SHE in the NM 
layer, which flows across the NM/FMI interface, depicted as 
the magenta arrow in Fig. |2}a). The spin orientation s of Js 
is oriented perpendicular to Js and Jq because of the SHE 
(cf. Eq.([TJ). At the FMI/NM interface, this spin current can 
exert a torque on and thus be absorbed by the ferromagnet if 
s is not parallel to the magnetization M of the ferromagnet. 
This absorption can only occur in a non-collinear configura- 
tion of M and s, since it only then is possible to transfer spin 
angular momentum from the spin current in the NM onto the 
magnetization of the FMI. In the Cartesian coordinate sys- 
tem defined by the unit vectors of the current direction (j), 
the surface normal (n), and the transverse direction (t = n x j, 
see Fig. |2c)) two different configurations are possible as il- 
lustrated by Fig. 12 a) and (b). In Fig. |2a) the magnetization 




FIG. 2. (Color online) Graphical illustration of the SMR exhibited by 
FMI/NM hybrid structures. A charge current Jq (blue arrow) flowing 
in the NM is converted via the SHE into a spin current Js (magenta 
arrow) flowing towards the FMI/NM interface. Due to the SHE the 
spin polarization s (violet arrow) is perpendicular to Jq and Js. At the 
FMI/NM interface the spin current is absorbed or reflected depend- 
ing on the relative orientation of the magnetization M in the FMI (red 
arrow) to s. The panels (a)-(b) show the two principle orientations of 
M with respect to the spin polarization s: If M is perpendicular to 
s it is possible to transfer angular momentum via the spin torque ef- 
fect, and the spin current in the NM gets absorbed by the FMI. For a 
collinear alignment between M and s, no angular momentum transfer 
is possible and the spin current is reflected at the interface, leading 
to a spin accumulation, (c) Definition of the coordinate system de- 
fined by j, t, and n in our YIG/Pt hybrid structures. In addition, the 
measurement scheme used for the determination /q, Vio^g, and Vtrans, 
with the Hall bar geometries w, I, and fpt is shown in the drawing. 



M||Jq is oriented along j and thus perpendicular to the spin 
polarization s of Js such that the spin current is absorbed. Fig- 
ure |2b) illustrates the case when M is oriented along t and 
thus collinear to s. In this case, the spin current cannot be 
absorbed by the FMI and is reflected at the interface. In the 
aforementioned Cartesian coordinate system, there is a third 
configuration where M is oriented along n. This case is analo- 
gous to Fig. 121 a), since M and s are again oriented perpendicu- 
lar to each other, resulting in an absorption of the spin current. 
In this simple picture we expect a higher resistance of the NM 
layer if the spin current is absorbed at the FMI/NM interface, 
and a lower resistance if not, such that Rm±s > ^mIIs- 

For a quantitative understanding of the SMR we need to 
calculate the spin diffusion process in a FMI/NM bilayer sys- 
tem and take into account the magnetization orientation de- 
pendence of the boundary conditions at the FMI/NM inter- 
face.~ These calculations are based on our first simplified 



theory presented in the supplements of our previous pubhca- 
tionr^ and show that the resuhing steady state charge current 
density Jq.finai is composed of two contributions: one along 
the initial current direction j and a second part along the trans- 
verse direction t. This leads to a magnetoresistance effect in 
the NM, which is sensitive to the magnetization direction of 
the FMI. The effect can be parameterized as a longitudinal 

Plong == Vlong Wfpt (^/q)"' and tranSVerse Ptrans = Vtransfpt (^q)"' 

resistivity (cf. Fig. |2c)) of the NM as a function of magneti- 
zation orientation m = M/Msat of the FMI^^i^ 



Plong = PO + Pl»I?, 
Ptrans = P2»In + P3»Jj»Jt- 



(2) 
(3) 
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According to the theoretical SMR model, p\ = — ps. |pi/po| 
is the SMR effect. p2 is a Hall-effect type resistivity. In 
addition to the conventional (ordinary) Hall effect of the Pt 
layer, p2 also contains a Hall-type SMR contribution.^" While 
Po, Pi, and p3 are expected to be independent of the exter- 
nal magnetic field magnitude, the Pt OHE coefficient con- 
tribution to P2 linearly changes^! with B, such that p2 = 
P2.SMR + P2.0HE(fi)- The magnetization parameters nij, nii and 
nia are projections of the magnetization orientation unit vec- 
tor m onto the coordinate system given by the unit vectors j, 
t, and n as illustrated in Fig.|2c). 

Following the calculations detailed in Ref. [Sfl the relative 
magnitude of the SMR effect is given by the ratios 
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Here, Anm is the spin diffusion length in the NM, Pnm the 
resistivity of the NM, Jnm the thickness of the NM, and G, 
the real part and G, the imaginary part of the spin mixing in- 
terface conductance. The relative magnitude of this effect, 
i.e., the SMR magnitude |pi/po|, is essentially determined by 
a|jj. Taking the Uterature value asH = 0.012 for NM = Pt 
(Ref. l45b. one thus expects a 10^^ relative resistance change 
for Pt. Moreover, the SMR effect will be large only if the 
thickness ^nm of the NM layer does not substantially exceed 
the spin diffusion length Anm in the NM. Last but not least, 
the SMR magnitude also will characteristically depend on the 
resistivity Pnm and the spin mixing interface conductance G, , 
at least as long as 2AnmPnmG,- < 1. 

It is important to compare Eqs.© and (|3]l to the anisotropic 
magnetoresistance (AMR) effect in conventional electrically 
conductive FMs^: 



Plong = Po + Apmj, 

Ptrans = P2'Wn + ApOTjmt. 



(6) 

(7) 



At first glance, AMR and SMR appear very similar In partic- 
ular, the SMR is not discernible from the AMR of a polycrys- 
talline ferromagnetic conductor if the magnetization resides 
in the ferromagnetic film plane (i.e., the j-t plane). This di- 
rectly follows from the unity relation 1 ~ nr +ml + m^, as 

we can rewrite m? into 1 — m^, for a magnetization oriented in 
the film plane («„ = 0). In contrast, if M has an out-of -plane 
component {nin ^ 0), we expect a different behavior: the mag- 
netization orientation dependence of the SMR vanishes for a 
rotation of the magnetization in the plane enclosed by j and n 
(oopt geometry, cf. Fig. |5jc)) since then m^ — 0. In contrast. 



the SMR depends on m for a rotation of m in the plane en- 
closed by t and n (oopj geometry in Fig. |5jb)) (mj — 0). For 
the conventional AMR of a polycrystalline FM the situation 
is reversed: no m dependence of the AMR in the t-n rotation 
plane (oopj geometry) and a clear m dependence in the j-n 
rotation plane (oopt geometry). The SMR thus differs quali- 
tatively from an AMR. 

In our experiments we measure the magnetoresistance as 
a function of the magnetization orientation m. More pre- 
cisely, in these so-called angle-dependent magnetoresistance 
(ADMR) experiments, one records piong and Pnans as a func- 
tion of the orientation of the externally applied magnetic field 
h = H/ |H|, while maintaining a fixed magnetic field magni- 
tude //meas = const.. In the very same fashion as in ADMR 
experiments in (Ga,Mn)As thin films /^'^'*''*^ we then use our 
SMR model equations to extract the po, pi, P2, ps parame- 
ters from the experimental data. This requires the knowledge 
of m for all field orientations h and field strengths //meas- To 
describe the magnetization orientation in the FMI we use the 
free enthalpy approach, with the expression 



GM{m) = -iJ.oH{h-m)' 



Bninl, 



(8) 



for the free enthalpy GmIhi), which takes an effective shape 
anisotropy contribution (Bn) and the Zeeman energy with the 
external magnetic field direction into account. We intention- 
ally neglect the (small) crystalline anisotropy of our YIG sam- 
ples, since we only consider experimental data taken at mag- 
netic fields much larger than the corresponding crystalline 
anisotropy fields. This allows to keep the number of free mod- 
eling parameters to a minimum. Using a numeric minimiza- 
tion algorithm, we determine the global minimum of the free 



enthalpy and assume that M always points along this orienta- 
tion. Please note that the free enthalpy in Eq.® is normalized 
to the saturation magnetization Ms. 



III. FABRICATION OF FMI/NM HYBRIDS 

If not explicitly indicated otherwise, the samples used in 
the experiments are YIG/Pt or YIG/NM/Pt heterostructures 
grown on (lll)-oriented gadolinium gallium garnet (GGG) 
or yttrium aluminium garnet (YAG) substrates (cf. Table H). 
The YIG films were epitaxially grown via laser-MBE from a 
stoichiometric polycrystalline YIG target, utilizing a KrF ex- 
cimer laser with a wavelength of 248 nm at a repetition rate of 
10 Hz4il^ A growth optimization yields the following opti- 
mum set of YIG growth parameters: 550°C substrate temper- 
ature during deposition, 2 J/cm^ energy density at the target, 
and an oxygen atmosphere of 25 /ibar. The NM layers were 
deposited in-situ, without breaking the vacuum, on top of the 
YIG film using electron beam evaporation, at room tempera- 
ture. 

After deposition, the structural and magnetic properties 
were investigated using high-resolution X-ray diffractometry 
(HRXRD) in a four circle diffractometer with monochromatic 
Cu Kai radiation, as well as superconducting quantum inter- 
ference device (SQUID) magnetometry. The results obtained 
for YIG/Pt hybrids grown on GGG substrates are compiled 
in Fig. |3ja)-(d). Figure [3ja) shows the sharp reflections of 
the GGG substrate and the laser-MBE grown YIG film in a 
20-(0-scan. The two low intensity reflections marked by as- 
terisks (*) are due to substrate impurities. The YIG reflec- 
tions, indicated by the arrow in Fig. |3jb) for the (444) plane, 
are only weakly discernible from the high intensity substrate 
reflection due to the small lattice mismatch of only 0.03% be- 
tween YIG and GGG. The high structural quality of the YIG 
layer is not influenced by the additional deposition of a Pt 
fibn as evident from Laue oscillations visible in the vicin- 
ity of the GGG (444) reflection in Fig. |3lb), which indicate 
a coherent, (1 1 l)-oriented growth of the YIG layer on GGG. 
From the position of the (444) reflection we calculate an inter- 
plane spacing of (^444 ~ 0.1787 nm from the Bragg equation. 
Assuming an undistorted cubic lattice structure this leads to 
a lattice constant of a = 1.238 nm, which is identical to the 
value flbuik of bulk YIG4^ Our result a = abuik is in contrast 
to published data on PLD grown YIG filmsi^"— In these pub- 
lications much larger lattice constants for the YIG layer have 
been reported, which result from a large rhombohedral distor- 
tion of the YIG lattice due to a deficiency of iron ions in the 
YIG structure i^2i^ The iron deficiency can be tuned by vary- 
ing the oxygen partial pressure during deposition.«24 From this 
we conclude, that our growth parameters allow the deposi- 
tion of highly stoichiometric YIG films from a stoichiometric 
polycrystalline target. 

To evaluate the mosaic spread of the laser-MBE grown YIG 
films we recorded X-ray rocking curves of the YIG (444) re- 
flection for our YIG film on GGG, one typical result is dis- 
played in Fig.[3jc). From a Gaussian fit (red line) to the data 
we extract a FWHM of 0.03°. This is an excellent value for 
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FIG. 3. (Color online) Structural and magnetic properties of the 
laser-MBE grown YIG/Pt hybrids, (a) 26-0) scan of a YIG (23 
nm)/Pt (7 nm) hybrid grown on a (1 1 l)-oriented GGG substrate. The 
two reflections marked with an asterisk (•) are imperfections present 
in the substrate prior to deposition, (b) The enlargement of the 
29 — CO scan in (a) reveals satellites due to Laue oscillations around 
the YIG (444) reflection, indicating a coherent growth, (c) The small 
full width at half maximum (FWHM) of 0.03° of the rocking curve 
of the YIG (444) reflection confirms the low mosaic spread of the 
YIG film and thus its excellent structural properties, (d) In-plane 
magnetization hysteresis curve of a 49 nm thick YIG layer on a GGG 
substrate determined from SQUID-magnetometiy at T = 300 K. The 
saturation magnetization Ms = 110 kA/m of the YIG layer is lower 
than the bulk value of Ms = 141.65 kA/m. The large noise is due to 
the large paramagnetic background signal from the GGG substrate 
(which was subtracted from the data in panel (d)). 



laser-MBE grown thin films and confirms the high structural 
quality of our samples. 

In the full range of the 20-o-scan (Fig.|3la)) we do not ob- 
serve any reflections which can be attributed to the Pt layer 
We therefore assume that our Pt thin films are polycrystalline 
and grow without a preferential texture on YIG. We evaluated 
the surface roughness of the laser-MBE grown YIG/Pt hybrids 
on GGG substrates via X-ray reflectometry and obtain from a 
simulation on average 0.7 nm for the surface roughness am- 
plitude of the YIG film and 0.8 nm for the Pt layer 

The YIG thin films deposited on diamagnetic, (111)- 
oriented YAG substrates (lattice mismatch; 3%) show similar 
structural properties. The 20-co-scans do not reveal any sec- 
ondary phases. The FWHM of the YIG(444) rocking curves 
is about 0.1° indicating a larger mosaic spread as compared 



to the films on lattice-matched GGG substrates. The X-ray 
reflectometry analysis gives on average a surface roughness 
amplitude of 0.8 nm for the YIG film and 0.9 nm for the Pt 
layer 

Taken together, our laser-MBE grown YIG/Pt heterostruc- 
tures exhibit excellent structural properties. Moreover, the 
structure of the YIG films is not influenced by the deposition 
of a Pt top layer 

The magnetic properties of our laser-MBE grown YIG films 
were analyzed using SQUID magnetometry. A hysteresis 
curve at r = 300 K is shown in Fig. [3td), with the external 
magnetic field H applied in the film plane. We subtracted 
the paramagnetic background signal of the GGG substrate and 
normalized the remaining magnetic moment to the volume of 
the 49 nm thick YIG layer to obtain the magnetization M of 
the YIG film. The M{H) hysteresis curve exhibits a low coer- 
cive field PLqHc = 3 mT and reaches a saturation magnetization 
Ms = 1 10 kA/m, which is approximately 80% of the reported 
bulk value Ms = 141.65 kA/m of YIG;^ We note that due 
to the large paramagnetic background signal the error for the 
determination of the saturation magnetization is large (at least 
10%) and thus alone might account for the measured differ- 
ence to the bulk value (cf. error bars in Fig. Eld)). Compared 
to the results reported by other groups , ^'^^^'^^•^^ we can con- 
firm that the saturation magnetization of YIG thin films on 
GGG substrates is close to the bulk value. However, the re- 
ported coercive fields in Refs. H^ H and E^ at T = 300 K 
are much lower (below 1 mT) than the one we observe in our 
YIG films (above 2 mT). On the other hand, the films de- 
scribed in Refs. IsT] and |53 exhibit coercive fields as large as 
15 mT. From these variations in the literature values for coer- 
cive fields and keeping in mind that H^ sensitively depends on 
the domain configuration, we conclude that the growth condi- 
tions greatly influence this quantity. 

For the YIG layers on (lll)-oriented, diamagnetic YAG 
substrates, we obtain a saturation magnetization of Ms = 
130 kA/m at r = 300 K, close to the bulk value of YIG. The 
extracted coercive field /io^c = 5 mT at T = 300 K for these 
samples is slightly larger than the values obtained for YIG 
films on GGG substrates. Nevertheless, these magnetic prop- 
erties confirm the state-of-the-art quality of our laser-MBE 
grown YIG thin films. 

An additional set of YIG/Pt heterostructures has been fab- 
ricated by sputtering a Pt layer on a commercially avail- 
able (FDK Corporation), 1.3 /im thick YIG layer deposited 
on (lll)-oriented GGG substrates via liquid phase epitaxy. 
The YIG film was grown under a PbO-B203 flux at around 
1200 K. For our liquid phase epitaxy prepared YIG/Pt het- 
erostructures we find in XRD studies a FWHM of the rock- 
ing curve for the YIG (444) reflection that is as narrow as the 
GGG (444) substrate reflection. Magnetization measurements 
on our liquid phase epitaxy YIG films yield Ms = 140 kA/m 
at r = 300 K close to the reported bulk value of YIGi^ 

The Fe304/Pt hybrid structure used in our experiments was 
also fabricated by laser-MBE and electron beam evaporation. 
First the magnetite layer was grown via laser-MBE on a (001)- 
oriented MgO substrate using a polycrystalline Fe304 target, 
with the following set of deposition parameters >S2j^ 320°C 



substrate temperature during deposition, 3.1 J/cm^ energy 
density at the target, and an argon atmosphere of 0.7 /ibar. 
After the magnetite growth, the Pt layer was deposited in- 
situ on top of the Fe304 film using electron beam evaporation, 
at room temperature. Structural characterization via HRXRD 
showed that no secondary phases were present in the samples. 
Satellites due to Laue oscillations around the Fe304 (004) re- 
flection and a narrow rocking curve of the (004) film reflection 
with a FWHM of 0.04° indicate the excellent structural qual- 
ity of these magnetite films .-^^^ SQUID magnetometry mea- 
surements^^i^ on these samples yield Ms = 450 kA/m close 
to the bulk value of 470 kA/m at T = 300 K (Ref.[60|). A Ver- 
wey transition at T = 117 K is visible indicating the perfect 
stoichiometry of the Fe304 thin film. 

The nickel ferrite (NiFe204) layer was grown via di- 
rect liquid injection chemical vapor deposition on (001)- 
oriented MgAl20/t^ at a substrate temperature of 700°C us- 
ing Ni(acac)2 and Fe(acac)3 (acac=acetylacetonate) as precur- 
sors. The Pt layer was sputtered ex-situ, after cleaning the 
surface of the nickel ferrite by Ar ion beam milling, at room 
temperature. The rocking curve of the NiFe204 (004) reflec- 
tion exhibited a FWHM of 0.2° measured via HRXRD. The 
magnetic properties were evaluated via an alternating gradi- 
ent magnetometer and yield a saturation magnetization^' of 
280 kA/m (bulk value'^^ ^^ ^ 399 kA/m). 



IV. ADMR EXPERIMENTS ON YIG/PT HYBRIDS 

For the magnetotransport experiments we patterned Hall 
bar mesa structures out of the plain, laser-MBE grown YIG/Pt 
films using photolithography and argon ion beam milling. The 
width w and length /hb of the Hall bar were 80 /im and 
1000 /im, respectively (Fig. |2tc)). We measure the magne- 
totransport at temperatures T between 5 K and 300 K in a 
liquid He magnet cryostat system in magnetic fields H of up 
to 7 T. The magnetoresistance was studied by applying a con- 
stant dc bias current density /q in the range of 0. 1 kA/cm^ to 
500 kA/cm^ along the Hall bar while recording the longitudi- 
nal Viong and transverse Vtrans voltage signals (Fig.[2|c)) of the 
Hall bar as a function of the external magnetic field magnitude 
H or orientation h. For the further evaluation we calculated 
the longitudinal piong = Viong(-/qO ' ^"d transverse resistiv- 
ity Ptians = Vtrans (-/qw)^', where / = 600 /im is the separation 
between the longitudinal voltage contacts and w = 80 /im the 
width of the Hall bar Within the range of charge current den- 
sities /q quoted above, the resistivities were independent of Jq. 
We note that a reference sample consisting of a single, laser- 
MBE grown YIG (30 nm) layer on a GGG substrate was found 
to be electrically insulating within our experimental limits. 
This gives a lower limit for the resistivity of our YIG layers 
Pyig > 6 X 10^ Hm, which is in agreement with the reported 
resistivity values for bulk YIG^ (Pyig = 6 x 10^ iim). 

In a first set of experiments we verified the existence of 
a magnetoresistance effect in our hybrids^^'22i2L^ by record- 
ing the resistance evolution as a function of the applied exter- 
nal magnetic field strength {R{H) measurements). FigureUa) 
and (b) show the R{H) results obtained aiT — 300 K from 
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FIG. 4. (Color online) R{H) measurements conducted on a laser- 
MBE grown YIG(54 nm)/Pt(7 nm) hybrid structure at T = 300 K 
grown on a (lll)-oriented GOG substrate, (a) Longitudinal resistiv- 
ity piong as a function of the applied in-plane magnetic field H with H 
oriented parallel (a = 0°) and perpendicular (a = 90°) to the current 
direction j. A clear resistance hysteresis behavior is visible in both 
orientations. The inset defines the magnetic field orientation angle a 
with respect to j. (The arrow indicates the positive angle direction.) 
(b) Hysteresis curve of the transverse resistivity for the external field 
applied at a = —45° and a = 45°. In both panels black symbols 
represent data taken during the //-field down-sweep, and red sym- 
bols represent data from the //-field up-sweep, as illustrated by the 
black and red arrows in the graphs. 



YIG(54 nm)/Pt(7 nm) bilayer grown on a ( HI )-oriented GGG 
substrate via laser-MBE. For the external //-field oriented par- 
allel to the current direction j (a = 0°) we observe character- 
istic hysteretic changes of pion^ around magnetic field values, 
that are larger than the coercive field values (cf. Fig.|3ld)). We 
attribute these differences to the Hall bar structuring process 
of the YIG/Pt hybrids. Upon changing the in-plane field orien- 
tation such that H is oriented perpendicular to j (a = 90°), the 
hysteretic behavior is inverted. Moreover, the resistivity val- 
ues for // > //c are clearly different for a = 0° and a = 90°. 
This dependence of the magnetoresistance on the external 
field orientation is expected from our theoretical model of the 
SMR (cf. Eq.©) and reflects the evolution of the magnetiza- 
tion orientation m = (mj,mt,»in) as a function of the external 
magnetic field magnitude in close analogy to the AMR effect. 
A similar behavior is expected for the transverse resistivity 
Ptrans(^)- From Eq.(|3]i the characteristic inversion should be 
most prominent for magnetic field orientations a = —45° and 
a = 45°. Our measurements in Fig. Ub) clearly confirm these 
expectations. Unfortunately, it is not straightforward to distin- 
guish directly between the polycrystalline AMR (Eqs.(|6]l, (jT)) 
and the SMR (Eqs.(|2|i, Q) effect from single R{H) measure- 



ments. Moreover, a quantitative evaluation of the MR effect 
from these type of measurements is difficult, since magnetic 
domain formation around He influences the observed MR sig- 
nal, while domain configurations are not taken into account in 
Eqs.dHi, ©orEqs.®, ©. 

ADMR measurements are not hampered by domain config- 
uration issues. In an ADMR measurement an external mag- 
netic field of constant magnitude //meas 3> //c is rotated in a 
given plane with respect to the sample, while piong and Po-ans 
are recorded as a function of this external magnetic field orien- 
tation. We here performed ADMR experiments in 3 rotation 
planes: ip, oopj, oopt. Herby, ip stands for a rotation of H 
around n (angle a), oopj for a rotation of H around j (angle 
/3), and oopt for a rotation of H around t (angle y). For the 
definitions of a, j3, and ysee the illustrations in Fig.|5la),(b), 
and (c), respectively, in which the field orientation h = H/H 
is also shown. Prior to the angular variation we initialize the 
magnetization by applying /io^init = 3 T along a = 0° (ip), 
j3 = Y= —90° (oopj, oopt). In our measurements we choose 
^meas to be at least one order of magnitude larger than the co- 
ercive field of YIG, such that the FMI is in a single domain 
state. 

We begin the evaluation of the SMR with the analysis of 
the ADMR results obtained from a laser-MBE grown YIG 
(54 nm)/Pt (7 nm) hybrid structure at T = 300 K. The experi- 
mental results are summarized in Fig.|5ja)-(c). For the ADMR 
experiments with an ip magnetic field rotation at /io^meas = 
1 T, /Xo//meas = 500 mT, and /io^meas = 100 mT, depicted 
in Fig. |5la), we observe an angular dependence with a pe- 
riod of 180° and a maximum in piong for h parallel (a = 0°) 
or antiparallel (oc = 180°) to j and a minimum for h parallel 
(a = —90°) or antiparallel (a = 90°) to t for every magnitude 
Miieas- Due to the small in-plane magnetic anisotropy of the 
YIG layer its magnetization is oriented parallel to the external 
magnetic field in good approximation for all a. This allows 
to directly use a to determine the magnetization orientation. 
Thus Piong is expected to follow a cos^ a dependence in ac- 
cordance to Eq.(|2]i, which is nicely reproduced by our ADMR 
data for all magnetic fields (lines in Fig.|5la)). In addition, the 
transverse resistance exhibits a cos a sin a dependence. We 
observe a maximum in Ptrans at a = —45° and a = 135°, 
while the minimum is located at a =45° and a = —135°. For 
both Piong and Puans the amplitude of the angular dependence 
is not influenced by the external magnetic field strength, and 
pi = — p3 as expected. 

In case of the oopj rotation plane (Fig. lUb)) we observe 
maxima in piong at /io//meas = 1 T and jLio//meas = 500 mT 
located at j3 = 0° (h || n) and j8 = 180° (h || -n). The 
minima in piong occur at j3 = -90° (h || t) and j3 = 90° 
(h II — t). At fioHmeas = 1 T the longitudinal resistivity again 
exhibits a cos^j3 dependence. Upon reducing //meas the shape 
anisotropy increasingly influences the orientation of the mag- 
netization in YIG, which results in deviations from the cos^ j3 
dependence's For the transverse resistivity we now observe 
a completely different angular dependence with a period of 
360°. Ptrans has a minimum at j3 = 0° (h || n) and a maxi- 
mum at j3 = 180° (h II — n). The ampHtude of this cos)3 de- 
pendence of Ptrans Strongly dependens on the external mag- 






FIG. 5. (Color online) ADMR results obtained from a laser-MBE grown YIG (54 nm)/Pt (7 nm) hybrid structure on a (lll)-oriented GGG 
substrate at T = 300 K. (a) Evolution of Piono (black squares) and Ptrans (red circles) as a function of h-orientation for an in-plane (ip) rotation 
of the external magnetic field at /io^meas = 1 T, ^o^meas = 500 mT, and /io^meas = 100 mT. The ip angle a is defined in the illustration above 
the data panels (the arrow indicates positive direction), (b) Dependence of piont; and ptrans on the oopj magnetic field angle jS at /Xo^meas = 1 T, 
Mo^meas = 500 mT, and /io^meas = 100 mT. The positive angle jS is illustrated in the drawing above the data plots, (c) Angular evolution of 
Plona and Ptrans for the oopt geometry at ^o^meas = 1 T, /io^meas = 700 mT, and ^o^meas = 100 mT. The oopt angle y is shown in the sketch 
above the experimental data. In panels (a), (b), (c), the green lines represent fits to the data using the SMR model defined by Eqs.©, O} with 
parameters of Table |l] while the dashed blue lines represent fits to the data using the AMR expression defined by Eqs.([6l(, (O. Interestingly, 
the angular dependence of Piong vanishes for the oopt geometry, which is not consistent with the AMR of a polycrystalline FM, but indeed is 
consistent with the behavior expected for the SMR effect. 



netic field strength (cf. Fig. |7]l. This is the signature of the 
ordinary Hall effect of the Pt layer The abrupt changes at 
Mo^meas = 100 mT, visible in piong and Ptrans. originate from 
the in-plane reorientation of the magnetization forced by the 
demagnetization energy, if the field is rotated near j3 = 0° or 
180°. 

Interestingly, the angular dependence of piong(7) vanishes 
for the oopt rotation plane (Fig. EJc)). This is in contrast to 
the expected angular dependence of an AMR of a polycrys- 
talline FM (cf. Eq.®), but is fully consistent with our SMR 
model (cf. Eq.©). Thus, we can exclude an AMR effect 
as the source of the observed MR. For Puans we again find 
a cos(7) angular dependence, which dominantly stems from 
the ordinary Hall effect in Pt. The abrupt changes in piong 
and Ptrans visible at /io^meas = 100 mT are explained by the 
abrupt in-plane reorientation of the magnetization, when the 
field orientation is rotated near 7 = 0° or y = 180°. 

Huang et al. also observed a MR in YIG/Pt bilayers^ and 
attributed it to an AMR effect originating from an induced 
magnetism (static magnetic proximity effect) in the Pt layer 



As discussed in the context of Fig. |5] our results suggest that 
systematic measurements as a function of M orientation in 
out-of-plane geometries allow to distinguish between SMR 
and AMR. The angular dependence we observe in our samples 
(Fig. Elc)) can be consistently explained by the SMR model 
presented in Sect. |II] while conventional AMR can be ruled 
out as an explanation for the observed magnetoresistance. 

For a quantitative analysis, we now employ the simulation 
technique successfully used for ADMR in (Ga,Mn)A o^^''*^'^^ 
and Heusler compounds^. First we choose a starting value 
for the anisotropy constant Bn and determine the magnetiza- 
tion direction m for every magnetic field orientation by nu- 
merically minimizing the free enthalpy given by Eq.(|8). We 
then calculate piong and Ptrans using m and a set of p,- param- 
eters via Eqs.(l2]i, (|3]i. An iterative optimization process of Bn 
and the p, parameters is carried out, until we achieve a sat- 
isfactory agreement {Y.X^ ^ 10^^) between experiment and 
simulation for all rotation planes and //meas magnitudes with 
a single set of £„ and Pi- For the simulation, all p, parameters 
except P2 were taken to be independent of the external mag- 
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netic field strength. The simulation curves obtained in this 
way are depicted as green lines in Fig.|5ja)-(c) and reproduce 
the experimental data very well. The respective p, parameters 
are summarized in Tableland we obtained £„ = 75 mT from 
the simulation. The p,- parameters extracted from our simula- 
tion show pi = — P3. This corroborates the prediction based 
on the qualitative description of the SMR effect. 

The MR of the laser-MBE grown YIG/Pt heterostructure 
only changes shghtly with temperature. At T = 5 K, pi/po 
is reduced by 34% as compared to T = 300 K. This decrease 
can be caused by a temperature dependence of the spin mixing 
conductance^^, the spin diffusion length, and the spin Hall an- 
gle. However, more experiments are required to further clarify 
this issue. Moreover, pi/po is independent of the current den- 
sity in the sample, such that thermal gradients generated via 
Joule heating can be excluded as a source of the MR effect. 

From the quantitative description of the SMR effect 
(Eq.llll)) it is evident that the SMR sensitively depends on the 
Pt layer thickness. We have thus used a set of samples with 
different Pt thicknesses and evaluated the SMR effect from 
ADMR experiments to extract the p, parameters in an itera- 
tion process as described above. The results of this procedure 
are compiled in Table|T] The ratio — pi/po is plotted as a func- 
tion of the Pt thickness in Fig. lUa). We clearly observe a 
maximum of the SMR in our samples (black symbols) at a Pt 
thickness of around 3 nmwitharatioof — pi/po = 1.6 x 10^^. 
A simulation of the experimental data using Eq.(|4|i allows 
us to extract the relevant material parameters from this film 
thickness dependence of the SMR. We however would like to 
emphasize that the parameter values thus obtained sensitively 
depend on ppt and the value used for G, . In the analysis, we 
explicitly took the film thickness dependence of ppt into ac- 
count by fitting Eq. dAll i to the Pt thickness dependence of po 
(see Fig. |9|a)) in the simulation. Note also that for values of 
fpt that are smaller than the surface roughness h, Eq. dAlb is 
no longer applicable, which puts a lower limit on the range of 
Pt thicknesses considered in the simulation. We note that /oo 
is the charge transport mean free path for an infinitely thick 
film (see Appendix [All : for a finite thickness the charge trans- 
port mean free path £ decreases with decreasing thickness.—. 
When surface roughness limits transport (fpt < /«,) our SMR 
model based on the spin diffusion equation is no longer ap- 
plicable. We can model such a regime by introducing a spin 
diffusion length or spin-flip relaxation time that depends on 
thickness of the metal layer^^"— A transport theory based 
on Boltzmann equations is better suited to treat the thick- 
ness dependence of the SMR in the presence of significant 
surface roughness scattering. However, such a more elabo- 
rate theoretical description is way beyond the scope of the 
present paper In particular, the solutions of a Boltzmann- 
type spin transport approach most Kkely will not be analyt- 
ical anymore. In contrast, the simple spin diffusion model 
we use here does yield analytical expressions, which allows 
to pinpoint the microscopic nature of the physical mecha- 
nisms. Realizing the limitations of our model, we focus on 
samples with fpt > 5 nm for the parameter optimization, since 
in this limit, the spin diffusion approach should yield a rea- 
sonable description of the spin transport processes. To cap- 
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no. 6. (Color online) Evolution of the SMR effect as a function 
of Pt thickness tpi. (a) fpt dependence of Pi/po determined from 
different laser-MBE grown YIG/Pt ([fpt] nm) samples (open sym- 
bols, squares for samples on GGG substrate, circles for samples on 
YAG substrates) from our fits to the ADMR data, (b) fpt dependence 
of the SMR signal for liquid phase epitaxy grown YIG and sputter 
deposited Pt heterostrucutures extracted from R{H) experiments^^, 
(c) Dependence of the MR on /pt determined by Huang et alM. for 
YIG/Pt hybrids. The lines in each panel represent a simulation based 
on Eq.lH for G,- = 4 x 10'"* Qr^m-- (red), G^ = 4 x lO'-^ Qr^mr'^ 
(green), and G;- = 4 x lO'^ fl^'m^ (blue). From this simulation we 
extract the spin Hall angle agn and the spin diffusion length Apt in 
Pt quoted in the respective panels (same color coding). 



ture the impact of G, we use three different typical values to 
show the influence of this parameter on the extracted asH and 
Apt. For Gr = 4 X lO'^ Qr^mr^ (cf. Jia et al.2i.) we obtain 
a satisfactory agreement between simulation and experiment 
(red line in Fig. |6ja)) for asH = 0.11, Apt = 1.5 nm. For a 
lower value Gr = 4 X lO'^ii^'m"', asH =0.25, Apt = 1.8 nm 
yield the best agreement between experiment and simula- 
tion (green line in Fig. HJa)). By choosing a higher value 
Gr = 4 X 10'^ Qr^mr^, we achieve for asH = 0.075, Apt = 
1 .2 nm the lowest deviation between experiment and simu- 
lation (blue line in Fig. |6ja)). However, the simulated SMR 
thickness dependence for Gr = 4 x lO'^ Qr^rcT^ diverges for 
small fpt, while the simulation for G,- = 4 x lO'^ H^'m^^ 
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sample 


T[K] 


po [lO-^flm] 


Pl/Po 


P3/P0 


P2(l T) [10-ilnm] 


GGG/YIG(54)/Pt(7) 
GGG/YIG(54)/Pt(7) 
GGG/YIG(54)/Pt(7) 


300 
30 

5 


406.5 
283.5 
208.7 


-5.9 xlO-'^ 
-4.7 X 10-4 
-3.9x10^4 


5.9x10-4 
4.7 X 10-4 
3.9 X 10-4 


-3.4 
-3.1 

-2.2 


GGG/YIG(61)/Pt(l.l) 


300 


1895.9 


-3.6 xlO-'^ 


3.6x10-4 


-2.5 


GGG/YIG(55)/Pt(1.2) 


300 


1857.3 


-3.9x10-4 


3.9 X 10-4 


-1.5 


GGG/YIG(57)/Pt(1.3) 


300 


1089.9 


-3.5x10-4 


3.5 X 10-4 


-18.0 


GGG/YIG(29)/Pt(1.8) 


300 


487.4 


-1.5x10-3 


1.5 X 10-3 


-2.5 


GGG/YIG(58)/Pt(2.2) 


300 


761.7 


-1.2x10-3 


1.2 X 10-3 


-3.5 


GGG/YIG(69)/Pt(2.7) 


300 


453.6 


-1.5x10-3 


1.5x10-3 


-4.8 


GGG/YIG(53)/Pt(2.5) 


300 


719.0 


-1.6x10-3 


1.6x10-3 


-5.8 


GGG/YIG(46)/Pt(3.5) 


300 


306.6 


-9.4 X 10-4 


9.4 X 10-4 


-3.0 


GGG/YIG(65)/Pt(6.6) 


300 


582.6 


-6.4 X 10-4 


6.4x10-4 


-7.8 


GGG/YIG(50)/Pt(7) 


300 


409.4 


-6.1x10-4 


6.1 X 10-4 


-3.4 


GGG/YIG(53)/Pt(8.5) 


300 


348.3 


-6.4 X 10-4 


6.4x10-4 


-3.5 


GGG/YIG(61)/Pt(ll.l) 


300 


334.5 


-4.4 X 10-4 


4.4 X 10-4 


-3.0 


GGG/YIG(52)/Pt(16.9) 


300 


339.2 


-3.2x10-4 


3.2 X 10-4 


-2.9 


GGG/YIG(55)/Pt(17.2) 


300 


331.7 


-3.0x10-4 


3.0x10-4 


-2.3 


YAG/YIG(45)/Pt(1.9) 


300 


1331.3 


-6.8 x 10-4 


6.8 X 10-4 


-3.5 


YAG/YIG(60)/Pt(2.5) 


300 


358.4 


-1.1x10-3 


1.1x10-3 


-2.8 


YAG/YIG(64)/Pt(3) 


300 


622.2 


-1.4x10-3 


1.4x10-3 


-3.2 


YAG/YIG(50)/Pt(3) 


300 


539.3 


-1.2x10-3 


1.2x10-3 


-4.0 


YAG/YIG(63)/Pt(6.5) 


300 


412.0 


-7.9x10-4 


7.9 X 10-4 


-3.4 


YAG/YIG(59)/Pt(6.8) 


300 


487.7 


-9.5 X 10-4 


9.5 X 10-4 


-3.5 


YAG/YIG(60)/Pt(9.7) 


300 


429.0 


-5.7 X 10-4 


5.7 X 10-4 


-3.9 


YAG/YIG(60)/Pt(12.8) 


300 


434.9 


-4.3 X 10-4 


4.3 X 10-4 


-3.5 


YAG/YIG(61)/Pt(19.5) 


300 


361.3 


-2.8 X 10-4 


2.8 X 10-4 


-2.9 


YAG/YIG(46)/Pt(21.7) 


300 


358.7 


-2.3 X 10-4 


2.3 X 10-4 


-2.7 


GGG/YIG(15)/Au(8)/Pt(7) 


300 


143.0 


-2.4 X 10-4 


2.4 X 10-4 


-11 


GGG/YIG(34)/Cu(9)/Pt(6) 


300 


111.0 


-0.9 X 10-4 


0.9 X 10-4 


-8 


MgO/Fe3O4(20)/Pt(7) 


5 


202.6 


-2.1x10-4 


2.1 X 10-4 


-4.2 


Mg AI2 04/NiFe204 (620)/Pt( 1 0) 


300 


242.0 


-2.7 X 10-4 


2.7 X 10-4 


-2.2 



TABLE I. The p, parameters obtained from fits to the experimental data for FMI/Pt hybrids and FMI/NM/Pt hybrids investigated in this work. 
The numbers in parentheses give the thickness of the respective layers in nm. All YIG based heterostructures in this table were fabricated via 
laser-MBE. The p,- parameters for the Fe304/Pt hybrid structure have been extracted at fl{)H = 1 T. 



and Gr = 4 X lO'^ ii-'m-^ nicely reproduces also the thick- 
ness dependence of the experimental data for small fpt. Since 
the value G, = 4 x lO'^ iJ-'m-" appears unreasonably high, 
when compared to the theoretical calculations by Jia et al.i^, 
we conclude that G^ = 4 x 10'4 Or'^vcT^ is most appropriate 
for our samples. The modeling for what we perceive to be 
reasonable G, values nicely reproduce the experimentally ob- 
served thickness dependence also in the limit of very thin Pt 
films (fpt < 5 nm). While the validity of the diffusive model 
is questionable in this limit, the Boltzmann corrections appear 
to be small. Nevertheless, a more elaborate theoretical evalu- 
ation is needed to clarify this point. 

In an independent set of experiments we investigated the 
Pt thickness dependence of the SMR in YIG/Pt heterostruc- 



tures fabricated via liquid phase epitaxy (YIG) and sputter- 
ing (Pt) on (lll)-oriented GGG substrates^^. These sam- 
ples were structured into Hall bars with w ~ 1000 /xm and 
^HB = 2200 /im. In this measurement series the longitudi- 
nal MR was determined from R{H) curves for each sample. 
The results of this analysis are shown in Fig. |6jb). For these 
samples we observe a maximum of the MR for a Pt thickness 
of around 5 nm with a MR ratio of 1.4 x 10-4. Following 
the same procedure as above, we included the thickness de- 
pendence of ppt via a fit of Eq. dAlb to the Pt thickness de- 
pendence of Po (see Fig. lUb)). For G^ = 4 x 10'4 Qr^rcC^ 
(red line in Fig. |6jb)) we obtain asH = 0.03, Apt = 2.5 nm, 
with satisfactory agreement between simulation and experi- 
ment. For a smaller Gr = 4 x lO'^ Qr^rcT^ (green line in 
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Fig-ISb)) the parameters change to asH = 0.06, Apt = 2.6 nm. 
A higher G,- = 4 x 10^^ Q.-^m-^ (blue line in Fig. I^b)) 
yields asH = 0.025, Apt = 2.2 nm. Again, the simulation for 
Gr = 4 X 10''*ii^'m^^ reproduces the experimental thickness 
dependence very well and is consistent with the theoretical 
calculations of Gr for YIG/noble metal interfaces. 

In the publication of Huang et alM the thickness depen- 
dence of the longitudinal MR in YIG/Pt heterostructures has 
also been investigated, the data taken from their publication 
are depicted in Fig.|6lc). Their MR thickness dependence has 
a maximum located also at a Pt thickness of 3 nm with a value 
of 5 X 10^^. A simulation of their data set with Eq.Q using 
ppt = 2.40 X 10^^ Hm (we use an average Pt sheet resistivity 
as no data on the thickness dependence is given in this pub- 
lication) yields Apt = 0.8 nm and asH = 0.11 for Gr = 4 x 
lO''* Qr^mT^, with an excellent agreement between data and 
simulation (red line in Fig.|6jc)). For Gr = 4 x lO'^ H^'m^^ 
the parameters in the simulation change to Apt =0.8 nm and 
asH = 0.34 (green Hne in Fig. I^c)). Gr = 4 x lO'^ Qr^mr^ 
(blue hne in Fig.|6tc)) gives Apt = 0.7 nm and a^n = 0.055. 
As we now used a thickness independent ppt for the simula- 
tion in Fig.|6tc), the simulation does not diverge for small fpt. 
While we anticipate that an inclusion of a thickness depen- 
dent ppt will affect the values for asH and Apt to some extent, 
the good agreement between the experimental data of Huang 
et al^ and our SMR simulation suggest that the SMR effect 
also plays a crucial role in their experiment. 

The spin diffusion lengths extracted from the three inde- 
pendent data sets in Fig. |6] compare reasonably well. The ex- 
tracted Apt depends on the assumed Gr and changes for our 3 
evaluated Gr values by 30%. While for a fixed Gr the relative 
error for Apt is less than 20%. Moreover, the spin diffusion 
length is comparable with the value of the charge transport 
mean free path in Pt (see Appendix|A]i. For fpt < 5 nm we find 
Apt > I, while Apt < £ is obtained for fpt > 5 nm. A Boltz- 
mann theory for the SMR at least at the level of the Fuchs- 
Sondheimer model will yield more insight into this problem. 
Also asH is strongly correlated with Gr and it is currently un- 
clear whether size effects enhance or reduce the SMR effect 
magnitude, leading to an effective, thickness dependent asH- 
The parameter values extracted from our analysis thus should 
be considered with this caveat in mind. We attribute the differ- 
ences in Apt and asH obtained for the different sets of samples 
to the different deposition techniques for the Pt layer Nev- 
ertheless, the values obtained from Fig. |6]deviate from those 
by Mosendz et al^ for Pt. Calculations^^ suggest that im- 
purities can substantially change the magnitude of the spin 
Hall angle, while they only slightly alter the spin diffusion 
length in Pt. In one of our recent publications we investigated 
the spin pumping effect in various conductive ferromagnet/Pt 
heterostructures^ and obtained a scaling relation, using the 
material constants quoted by Mosendz et alJ^^ for Pt. We 
would like to note that the scaling relation also holds for the 
Pt parameter values we extracted from the SMR experiments 
in Fig. |6] Using asH = 0.11 and Apt = 1.5 nm we obtain a 
spin mixing conductance in conductive ferromagnet/Pt bilay- 
ers of Gr = 4 X lO'"* Qr^rcT^, corroborating the analysis of 
the SMR experiments above. 



We now address the Hall effect in our laser-MBE grown 
YIG/Pt hybrids in more detail. The field dependence of p2 
for T = 300 K and T = 30 K extracted from the simulation 
of our ADMR data is shown in Fig. [T] for the very same 
sample as shown in Fig. [5] (YIG (54 nm)/Pt (7 nm)). We 
here plot P2 versus the magnetic inductance lio{H +Myig) 
to take into account the additional magnetic field in the Pt 
layer due to the magnetization of the YIG. The analysis is 
only conducted for /Xo//meas > 0.25 T as only then the mag- 
netization of the YIG layer is saturated and oriented along 
the external magnetic field direction. For T = 300 K we 
use Myig = 110 kA/m determined from SQUID magnetom- 
etry (Fig. lUd)). For T = 30 K we extrapolated the satura- 
tion magnetization from 300 K using the temperature depen- 
dence of bulk YIG^ and obtain Myia = 150 kA/m. A di- 
rect determination of the saturation magnetization via SQUID 
magnetometry at this temperature is not possible because of 
the paramagnetism of the GGG substrate. At both temper- 
atures p2 increases linearly with increasing magnetic field 
(/Zo(// +Myig)), as one would expect for an OHE.-** A lin- 
ear fit to the data to extract the Hall coefficient of our Pt 
yields rnaii ~ —2.5 x 10^'^ m^/C for both temperatures. This 
value is close to —2.1 x 10^" m^/C reported in Ref. l73lfor 
evaporated Pt films. Moreover, we find a non vanishing ab- 
scissa Pahe = (-0.61 ±0.1) X 10"" i2m at r = 300 K and 
Pahe = (-0.41 ±0.1) X 10-" i2m at r = 30 K. This effect 
can be understood within our theoretical modeP" of the SMR 
effect, which predicts an anomalous Hall-like SMR contribu- 
tion due to the imaginary part of the spin mixing conductance 
(cf. Eq®). From our experiment we obtain for fpt = 7 nm a ra- 
tio Pahe/Po = 1 .5 X 10^^ at r = 300 K. Using the parameters 
Gr = 4 X lO'"* i2-im-2, asH = 0. 1 1, Ar = 1 .5 nm and Eq.© 
we extract G,- = 1 . 1 x lO'^ ii^ 'm^^ from the experiment. This 
gives Gi/Gr = 0.03 which nicely agrees with theoretical cal- 
culations^ {Gi/Gr ~ 1/20). The quantitative agreement be- 
tween theory and experiment for both magnetoresistance and 
Hall-type measurements confirms the existence of the SMR 
effect. Note that an induced ferromagnetism at the YIG/Pt in- 
terfac o^^'^^ — might also contribute an AHE signal. However, 
our quantitative analysis of the Hall data, together with the 
MR data with out-of-plane magnetization have little room for 
such a static, magnetic proximity induced MR effect. 

While the magnetoresistive behavior observed in our 
YIG/Pt samples is not consistent with an induced magneti- 
zation (static magnetic proximity effect) in the Pt close to the 
YIG/Pt interface, an experiment to rule out magnetic proxim- 
ity as the origin of the SMR is desirable. In a recent publica- 
tion^ we employed XMCD measurements to detect element- 
specific magnetic moments at the Pt Ls-edge in the very same 
laser-MBE grown samples also used for SMR measurements 
presented here. The XMCD measurements showed that the 
induced magnetic proximity moment in Pt - if present at all - 
is small in our samples. However, the available XMCD data 
do not allow to exclude a static proximity effect in Pt, only an 
upper limit can be put on the induced Pt moment. Therefore, 
another set of experiments addressing the magnetic proximity 
effect in Pt is desirable. 

Following the same line of argument as in our previous 
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FIG. 7. (Color online) Magnetic inductance dependence of P2 for a 
laser-MBE grown YIG (54 nm)/Pt (7 nm) hybrid structure extracted 
from the fit to the ADMR data at T = 300 K (black open squares) 
and r = 30 K (blue closed circles). The black and blue lines repre- 
sent linear fits to determine the Hall coefficients rnaii and the AHE 
contributions Pahe for T = 300 K and T = 30 K, respectively. 



publication on the SMR,— we also investigated the impact of 
a NM layer between YIG and Pt on the SMR. In Fig. [ga)- 
(d) we show the ADMR signals of a laser-MBE grown YIG 
(45 nm)/Au (7 nm)/Pt (7 nm) and a laser-MBE grown YIG 
(45 nm)/Cu (9 nm)/Pt (7 nm) heterostructure. We note that 
the discrepancy between the green simulation curve and the 
experimental Ptrans data in the ip rotation for both samples 
arises from a small out-of-plane external field component 
caused by a slightly tilted rotation plane (2°). This experi- 
mental misalignment leads to a superimposed cos a depen- 
dence. As the OHE and AHE contribution p2 is compara- 
ble to p3 in these samples, a slight tilting in the experiment 
from the ideal ip rotation plane already leads to a considerable 
deviation between experiment and simulation. The ADMR 
signal (Fig.[8la)-(d)) from the laser-MBE grown YIG/NM/Pt 
samples is qualitatively identical to the ADMR from YIG/Pt 
(Fig.|5]l. ADMR experiments in YIG (45 nm)/Au (7 nm) and 
YIG (45 nm)/Cu (7 nm) samples without Pt top layers did 
not reveal any variation of the resistance with H oriented in- 
plane within our experimental resolution. From these exper- 
iments, we obtain an upper limit of — pi/po < 5 x 10^^ for 
YIG/NM samples without Pt (data not shown here). As men- 
tioned in Sec.inithe maximum — pi/po is strictly smaller than 
a^^. Assuming asH = 0.0035 for Au (Ref. 1451). one estimates 
— Pi/po < 1.2 X 10^^, which presently is below our experi- 
mental resolution. Thus, we are currently unable to resolve 
the SMR in FMI/Au hybrids. These results demonstrate that 
the MR effect observed in FMI/Pt and FMI/NM/Pt samples in- 
deed is related to spin currents. In particular, a static magnetic 
proximity effect, i.e., an induced magnetic moment in the Pt 
layer, can be ruled out based on these experiments, since a 
magnetic proximity effect can not persist over 8 to 9 nm of 
Cu or Aui2^ For the SMR, only the conversion of a charge 
current into a spin current and back into a charge current via 
the SHE/ISHE in the Pt and the transport of spin current from 
the Pt to the YIG are necessary. Thus, the SMR persists even 
when a NM layer is inserted, provided that the spin current 



can propagate through this NM layer The extracted quanti- 
tative data for all heterostructures from the ADMR simula- 
tion are summarized in Table |I] po decreases with increasing 
total NM/Pt layer thickness, since in these double NM struc- 
tures the two parallel conducting layers both contribute to the 
total resistivity and since the Au and Cu layer have higher 
conductivity. The two-layer parallel conductance is also ev- 
idenced by the increase in p2(l T) compared to the YIG/Pt 
hybrid structure, due to the larger absolute Hall constant of 
Cu and Au.^SiSS, The ratio Pi/po decreases by a factor of 3 for 
the YIG/Au/Pt heterostructure compared to the YIG/Pt ref- 
erence sample. For the YIG/Cu/Pt heterostructure the ratio 
even decreases by a factor of 7. This decrease in the SMR 
effect can be rationalized in terms of the exponential decay 
of the spin current determined by the spin diffusion length in 
the NM and the parallel conduction channel which NM repre- 
sents. A quantitative theoretical description would be highly 
desirable, since it would possibly allow to use SMR experi- 
ments in YIG/NM/Pt hybrids with varying NM thickness to 
extract the spin diffusion length in more complex structures 
using simple ADMR experiments in the future. 

Finally, we show that the SMR effect is not limited to 
FMI/NM heterostructures based on the ferromagnetic insu- 
lator YIG. According to our model any FMI/NM bilayer in 
which the resistance of the FMI is several orders of magnitude 
larger than the resistance of the NM layer (such that the resis- 
tance of the NM dominates) should exhibit a SMR. To verify 
this conjecture we also investigated samples based on other 
ferromagnetic insulators and semiconductors. One sample 
consists of a 20 nm thick, (OOl)-oriented magnetite (Fe304) 
layer on a (OOl)-oiiented MgO substrate covered in-situ by an 
electron beam evaporated, 7 nm thin, Pt film. We note that 
the resistivity— of Fe304 is of the order of 1 x 10^ niim at 
T = 300 K, which is two orders of magnitude larger than the 
sheet resistivity of the Pt layer Therefore, we performed the 
ADMR measurements at T = 5 K below the Verwey transi- 
tion, where Fe304 becomes semi-insulating^ ~ such that the 
conductance is dominated by the Pt layer (Fig. [8je,f)). We 
note that the observed SMR ratio and the respective p, pa- 
rameters for our Fe304/Pt heterostructure depend on the ex- 
ternally applied magnetic field, which is caused by the high 
saturation field of Fe304 due to antiphase boundaries in the 
filmr^Sli Another independent sample consists of a semicon- 
ducting, 620 nm thick nickel ferrite (NiFe204) layer on a 
(OOl)-oriented MgAl204 substrate with an ex-situ, sputter de- 
posited 10 nm thick Pt layer (Fig. |8lg,h)). The resistivity of 
these NiFe204 films is 1x10^ niim at T = 300 K with a 
bandgap of 0.5 eV. We then performed ADMR experiments in 
ip, oopj, and oopt geometry as described above. In Fig. ISje)- 
(h) we show the corresponding results for the ip and oopj ge- 
ometry. The parameters of the fits extracted from the exper- 
iments at ^.qH = 1 T are summarized in Table J] In all cases 
the fits reproduce very well the angular evolution of the data. 
The MR observed in experiment thus is consistent with SMR, 
but inconsistent with AMR. Taken together, we thus have ob- 
served the SMR effect in all FMI/Pt and FMI/NM/Pt samples 
studied. 
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FIG. 8. (Color online) Magnetoresistance data (symbols) and corresponding SMR fits to Eqs.©, (|3]l (green lines) from laser-MBE grown 
YIG (15 nm)/Au (8 nm)/Pt (7 nm) (panels (a) and (b)), and laser-MBE grown YIG (34 nm)/Cu (9 nm)/Pt (7 nm) (panels (c) and (d)), Fe304 
(20 nm)/Pt (7 nm) at T = 5 K and yij^H = 1 T (panels (e) and (f)), NiFe204 (620 nm)/Pt (10 nm) (panels (g) and (h)) at T = 300 K and 
\i{)H = 1 T. The oopt rotations (not shown here) exhibit no angular dependence in Piong for all 4 samples. The introduction of a second NM 
layer between YIG and Pt influences only the magnitude of the observed SMR. This clearly indicates that the SMR arises due to the spin 
current absorption at the YIG/NM interface and is not due to the MR behavior of an induced magnetic polarization in the Pt layer. The SMR 
signature also appears in other FMI/Pt hybrid structures and thus confirms the universality of the SMR effect. 



V. CONCLUSIONS 



In summary, 'we have quantitatively investigated the SMR 
effect in FMI/Pt and FMI/NM/Pt heterostmctures. The SMR 
effect is based on the conversion of a charge to a spin cur- 
rent via the spin Hall effect, and back to a charge current 
via the inverse spin Hall effect in the NM (Pt) layer The 
SMR effect characteristically depends on the absorption of 



the spin current at the FMI/Pt or FMI/NM interface, 'which 
in turn can be tuned via the orientation of the magnetiza- 
tion of the FMI with respect to the spin polarization of the 
spin current. Thus, the SMR effect enables a remote sensing 
of the magnetization direction in the FMI by simply measur- 
ing the resistance of the adjacent NM layer We have shown 
that the signature of SMR is qualitatively different from con- 
ventional AMR, in particular when the magnetization has a 
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component perpendicular to the Pt film plane. Magnetotrans- 
port measurements as a function of the magnetization orienta- 
tion, rotating the magnetization from within the sample plane 
to a perpendicular-to-plane orientation, thus allow to disen- 
tangle AMR and SMR. We have observed SMR in YIG/Pt 
hybrids, as well as in YIG/Au/Pt, YIG/Cu/Pt, Fe304/Pt and 
NiFe204/Pt heterostructures. These results confirm the SMR 
thus as a universal and robust effect, that is not limited to cer- 
tain material combinations. This allows to use the SMR in 
a wide variety of material combinations. Moreover, the ex- 
cellent quantitative agreement between theory and experiment 
clearly show that a static magnetic proximity effect is not the 
origin of the observed MR in FMl/NM hybrids. We demon- 
strated that our quantitative model for the SMR effect allows 
us to extract the spin Hall angle asH =0.11 ±0.08 and the 
spin diffusion length Anm ~ (1 .5 ± 0.5) nm of Pt from SMR 
experiments, by examining the thickness dependence of the 
SMR effect. This enables to use the SMR for the determina- 
tion of the spin diffusion length and spin Hall angle in various 
NMs. The SMR is a novel, simple to measure magnetoresis- 
tance effect, which paves the way for new spin current related 
experiments. Last but not least, our SMR theory predicts a 
small AHE-like SMR contribution to the transverse resistivity 
due to the imaginary part of the spin mixing conductance. In 
our experiments we find Pahe/Po = 1 -5 x 10^^ atT ~ 300 K 
for fpt = 7 nm. 



ing at the interfaces, for our analysis we strictly use the dif- 
fusive limit {p = 0). In FigUJa) we plotted the po(fpt) for 
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Appendix A: Thickness dependence of Pt resistivity 

For the simulation of the thickness dependence of the SMR 
ratio it is necessary to also include the thickness dependence 
of ppt. The thickness dependence of the resistivity for thin 
films has already been extensively studied in experiment and 
theory.-SSrJ^ For our analysis we use the expanded version 
of the Fuchs-Sondheimer theory , ^^'^^ that includes a surface 
roughness amplitude described by the parameter h. In the 
limit that fpt > h, the thickness dependence of ppt can be writ- 
ten asi^ 



PPt(fPt)=Poo 1 



8itpt-h) 



iU^-p)) 



(Al) 



where Poo is the resistivity, and [„ the mean free path for an 
infinitely thick film. The parameter p describes the scatter- 



FIG. 9. (Color online) Evolution of the Pt resistivity as a function of 
Pt thickness, (a) Platinum thickness dependence of po determined 
from different laser-MBE grown YIG/Pt ([tpi] nm) samples (open 
symbols) from our fits to the ADMR data (cf. Fig.|6ja)). (b) Pt thick- 
ness dependence of the resistivity for liquid phase epitaxy grown 
YIG and sputter deposited Pt heterostrucutures (cf. Fig. [Sfb)). The 
red line in each graph is a fit to the experimental data via Eq. dAlt . 

different laser-MBE grown YIG/Pt ( [fpt] nm) samples (cf. Ta- 
ble m. As expected the resistivity increases with decreasing 
film thickness. By fitting the experimental data by Eq. dAll i we 
obtain the following set of parameter: paa = (3.35 ±0.63) x 
10-^ iim, loo = (4.1 ±2.0) nm, h = (0.8 ±0.1) nm. Some 
experimental values deviate from the fitted theoretical curve 
as the assumption of a thickness independent surface rough- 
ness might be to simplistic. Nevertheless, the fitted curve re- 
produces the trend from our experiment and thus should be 
sufficient enough for the simulation of the SMR ratio. More- 
over, the value for h agrees nicely with the average surface 
roughness amplitude extracted from X-ray reflectometry mea- 
surements (cf. Section HiH) . The results we find for the LPE 
grown YIG/Pt hybrids are depicted in Figj9jb). Here a fit to 
the experimental data gives: poo = (7.00 ±0.65) x lO^^iim, 
/oo = (2.1 ± 1.1) nm, h = (1.8 ±0.2) nm. These two sets of pa- 
rameter enable us to insert Eq. dAlt into Eq.(|4) to describe the 
thickness dependence of ppt in our simulation in Fig.|6ta) and 
(b). Taking into account the different deposition techniques, 
these values compare reasonably well to published values by 
Fischer et al^ (poo = 160 niim, /oo = 10 nm) for Pt thin films 
evaporated onto glass substrates. 
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